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Abstract 

We calculate finite volume effects on the pion masses and decay constant in twisted mass lattice 
QCD (tmLQCD) at finite lattice spacing. We show that the lighter neutral pion in tmLQCD gives 
rise to finite volume effects that are exponentially enhanced when compared to those arising from 
the heavier charged pions. We demonstrate that the recent two flavour twisted mass lattice data 
can be better fitted when twisted mass effects in finite volume corrections are taken into account. 
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I. INTRODUCTION 



For large enough lattices, finite volume effects (FVEs) are only sensitive to the long dis- 
tance physics of the underlying theory - for lattice QCD simulations they are well described 



by chiral perturbation theory (ChPT) and can be analyzed systematically in a chiral 
expansion. In a series of papers |2"^, FVEs on pseudoscalar meson masses and decay con- 
stants in the continuum limit have been studied with the help of resummed Liischer formulae 
used in conjunction with ChPT. In the original asymptotic formula by Liischer jsl, the expo- 
nentially dominating FVEs for the mass of a particle P is expressed as an integral over the 
forward scattering amplitude of the particle P off the lightest particle in the spectrum (pions 
for QCD). Since the scattering amplitude is needed only at low energy, it can be evaluated 
in ChPT. In Ref. {4], it has been shown that by resumming a series of such integrals (over 
the same scattering amplitude but with kernels increasingly exponentially suppressed), one 
can improve the reliability of such formulae in describing FVEs. 

As shown in Ref. j^, if one inserts the tree- level scattering amplitude in the resummed 
Liischer formula, one obtains exactly the one- loop ChPT calculation of the FVEs. In Ref jo], 
the pion mass has been calculated in ChPT to two loops and compared with the result 
obtained from the resummed Liischer formula using the one-loop ChPT representation of the 
scattering amplitude. The difference is found to be very small, and one has thus confidence 
in the validity of the resummed Liischer formulae in predicting the size of FVEs. This 
allows for a much simpler evaluation of the main effects at the two-loop level and beyond, 
and for one to check the convergence of the chiral expansion for FVEs. Nevertheless, it is 
still essential to test these predictions against actual data. 

Recently, the European Twisted Mass Collaboration (ETMC) has provided one such test 
for the pion mass and decay constant 0, Isl, and the outcome was not entirely positive. 
While FVEs on the decay constant are rather well described by the resummed formula at 
the next-to-leading order (NLO), the same does not hold for the pion mass. For the pion 
mass, FVEs are larger than what is predicted from the NLO and NNLO resummed formulae 

n 

(see Table 2 in Ref. [7j]), and the discrepancy appears to be much larger than the size of the 
error estimated in Ref. Q], which calls for an explanation. 

Before discussing possible sources of such a deviation, we stress that the ETMC data are 
on the borderline of applicability of ChPT for calculating finite volume corrections. It has 
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been argued in Ref. ^ that for ChPT to be applicable, a box size of L > 2 fm is necessary. 
The comparison made in Table 2 of Ref. (tI uses the largest volumes as reference points, 
and FVEs are measured with respect to these. The data sets used to study FVEs thus have 
L ~ 2.0 fm (where ChPT may still marginally work) or less (where ChPT is not expected 
to be valid). This may be a possible reason for the discrepancy. 

In this paper we study another possible source for deviation from the continuum formulae, 
i.e. discretization effects. As is well known, at finite lattice spacing, isospin- and parity- 
breaking effects are sizeable in twisted mass lattice QCD (tmLQCD). In particular, the 
neutral pion mass becomes smaller than the charged ones ^, and parity-breaking interactions 
among pions become possible. Both of these effects are well described in the framework of 
tmChPT jo-JJ], and they have nontrivial influences on the finite volume corrections as 
they generate exponentially enhanced FVEs. Indeed, if one takes the continuum limit first 
without explicitly accounting for these effects and make finite volume corrections using the 
continuum formulae, one will not be able to fully disentangle in the final results FVEs from 
discretization effects.^ 

The rest of the paper is organized as follows. In the next section, we briefly describe 
the formalism of tmChPT, highlighting aspects that are relevant for our calculations. In 
Sec. mil we explain why and which discretization effects are exponentially enhanced at finite 
volume with the help of LO asymptotic formulae. In Sec. llVt we give resummed asymptotic 
formulae at NLO applicable to the ETMC data, which we use to perform a new analysis on 
the ETMC data. Sec. |V] contains relevant results and discussion of our analysis, and Sec. NH 
our conclusion. We concentrate in this paper on the charged pions where high quality lattice 
data are available, but also provide formulae for FVE for the neutral pion mass for future 
use. 



^ For the range of masses and lattice spacings of the ETMC simulations, the neutral pion is typically about 
15-20% lighter than the charged ones. 

^ The size of the twisted mass discretization effects can only be determined after a detailed analysis. Nev- 
ertheless, one can argue that if these were responsible for the discrepancy between the FVEs observed by 
ETMC and those calculated in the continuum, they are at the percent level. 
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II. TWISTED MASS CHIRAL EFFECTIVE THEORY 



Consider tmLQCD with a degenerate doublet of quarks [l3|. The low energy, long dis- 
tance dynamics of the underlying lattice theory can be describe d by an effective continuum 
chiral theory constructed using the two-step procedure of Ref. [ij], whereby effects of dis- 
cretization errors are systematically incorporated in a joint expansion of the lattice spacing, 
a, the quark mass, m, and the twisted mass, fi. This was carried out to NLO in Ref. 
and the resulting tmChPT studied in detail in Ref. 12|. Using the power counting scheme, 
m ~ yU ~ ~ (2Aq(-.j-), the twisted mass effective chiral Lagrangian reads 12|: 



-^j{D,UD,Ur-^j{D,UDM'? 
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(x'^u + u^xf 



- w{x'^u + u^x') {A^u + u^A) - w'{A^u + u^Af 

+ WioiD^A^D^U + D^U^D^A) + W{D^UD^U^){A^U + A) 



(1) 



where we have displayed only parts relevant for our discussion and our study of FVEs below. 
Here, F is normalized so that = 92.4 MeV, U is the usual SU(2) matrix-valued field, 
(. . .) denotes the trace, ii the usual Gasser-Leutwyler Low Energy Constants (LECs), and 
the rest of coefficients LECs arising from discretization errors. The quantities x' ^md A are 
spurious for quark masses and discretization errors that are set to constant values at the 
end of the analysis: 



x' 2Bo{m + ir3/i) + 2Wo a = m + a + ir^fi , A 2Wo a 



(2) 



where is normalized so that r| 



1, Bq and Wq are unknown dimensionful constants, and 



we have defined the quantities m, u, and a. 

n 

As explained in Ref. [Ij], the Wiq term in the effective chiral Lagrangian is redundant, 
and can be transformed away into a combination of W, W and H' terms. However, the same 
redundancy can be used to eliminate the W term in favour of the Wio term instead, which 
has the advantage of simplifying the Feynman rules when studying FVEs in tmChPT, as 



we see below. In doing so, results derived in Ref. 



12l | would remain the same except with 



all W terms removed. In particular, the expansion about the NLO vacuum with external 
fields set to zero now reads: 



= 1^- ■ d,if +^{M' + AM') - 4^,aW + - (vr^)^ 



M' 



+ 



(3) 



- -p - ■ 2 ■ ' " ' 2 F2 2 F ' ' 24F2 

where vr^ = vf • vf, and we have used the convenient parametrization U = a + m ■ t/F, with 



o" = a/1 — TT^/F^. The quantities s and c denote the sine and cosine of a nonperturbatively 
defined vacuum (twist) angle u, and e is the shift of the vacuum angle at NLO tmChPT 
from that at LO: 

16 

e = uJm-uo = -j-2^So {W + 2W'acolM') . (4) 

Note that uj differs from both loq and Um by 0{a) so that at NLO accuracy, s and c could 
equally well be Sq (sm) and Cq (cm)- The mass parameters are given by 



M' 



'I = V(m + a)2 + /i2, AM' = — [M'^fg + lQ{acM'W + a^c^VT 



(5) 



and the charged and neutral pion masses at NLO read: 

M' ^ 32M' 



Mi 



M' 



327r2F2 F2 
32 



ac(M'iy + acW 
128 



with 



M^o = Ml± - -^ah^W = Ml± -a^K, K= —s^W^W . 



+ 2 log — — = log — ^ 

^ ^ M. . ^M2 



(6) 
(7) 

(8) 



the standard scale-independent LECs of SU{2) ChPT [15]. We have defined the dimensionful 
quantities Aj for use in our numerical analysis below. The pion decay constant at NLO reads 



J. 77" -L -jT- 



M' - AM' 



(9) 



167r2F2 

Note that at maximal twist, u = 7r/2, the charged pion masses and pion decay constant 
above take their continuum form with 0{a) effects automatically removed. At NLO, the 
effects of twisting reside only in the O(a^) charged-neutral pion mass splitting, which is 
maximal at maximal twist. 
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III. EXPONENTIALLY ENHANCED TWISTED MASS DISCRETIZATION EF- 
FECTS AT FINITE VOLUME 



Consider Liischer's formula for the pion mass ^: 

oo 

M^{L) - 3 



j dyT{iy)e-^'^- + . . . , K = M^L, (10) 



where Tiiy) is the forward vrvr-scattering amphtude evaluated for a purely imaginary ar- 
gument (given in units of pion mass), and y is real and dimensionless. This formula gives 
the exponentially dominant contribution under the assumption that the pion is the lightest 
particle in the spectrum. In the presence of a splitting between the charged and the neutral 
pions, as is the case in tmLQCD, this formula is modified to 



_ M^,(L)-M 
tii = — 



3 

-^^j^yr^iM^^^'^-, (u) 



167r2A 

where J^ij{iy) is now the forward scattering amplitude of pions with isospin index i off pions 
with isospin index j. Note that y is still dimensionless but now normalized to M^j, and we 
use the shorthand Xj for X^j. 

With its lighter mass, the neutral pion contribution is (exponentially) dominant with 
respect to the heavier charged pions. The weights of these contributions are given by the 
forward scattering amplitudes, and at LO in tmChPT they read: 

M' 

•^ll = --Fl2 = --^13 = ^, (12) 

where M' and F are the LO pion mass-squared and decay constant. This shows that since 
pions with index 1 and 2 are degenerate, their contributions cancel. The relative finite- 
volume shift for the charged pion mass at LO in tmChPT is then given by 

where Ki denotes the modified Bessel function. As an illustration of the importance of 
this discretization effect, we evaluate R± for the ETMC Bi ensemble in Ref. j?], for which 



For clarity in demonstrating the idea here, we do not write out the resummed version of this formula, but 
will defer to the next section where we give our final formulae in full. 



6 



M^± = 0.33 GeV and M^o = 0.27 GeV: 

R± = 0.26% (Bl ensemble) . (14) 

However, if we set M^o = M^± = 0.33 GeV, i.e. turning off the twisted mass effects, we get 

= 0.15% (Bl ensemble) . (15) 

We note that if one here resums the whole series of exponentially subdominant LO contribu- 
tions (corresponding to a LO ChPT evaluation of FVEs), 0.15% goes up to 0.62%; using the 



NNLO resummed asymptotic formula of Ref. j4| gives 1%, whereas the ETMC measurement 
is 1.8(5)%. 

There is a second discretization effect arising from the parity-violating cubic interactions 
(see Eq. ([3])), which could be potentially significant as well. The contributions to FVEs 
due to these in the pion mass formally come in at higher order: they are NNLO FVEs, or 
0{p^) corrections to the pion mass.^ But because of the different topology of the relevant 
loop diagrams, these contributions are exponentially enhanced with respect to the tadpole 
diagrams that contain virtual neutral pions. For these diagrams, the dominating exponential 
behavior goes as e'^"^^^^^, where w = M^o/(2M^±). From the parameters extracted from 
the Bl ensemble, the enhancement of such a contribution with respect to the standard e~^='= 
behavior is more than 200%, i.e. e^^°''^^^^ /e^^^ ~ 2.4, which motivates us to investigate 
these effects in our analysis even when they are formally of higher order. Now these con- 
tributions to the FVEs are proportional to e^, which involves unknown LECs W and W 
that have to be determined if the size of these contributions are to be predicted. If one 
works at maximal twist, only W is required, which one can determine e.g. from the ratio 
of density matrix elements jl2|.^ We remark here that if the exponentially enhanced parity- 
violating contributions are truly large or that FVEs can be measured precise enough, one 
can turn it around and use instead FVEs to get a measure of W; one could even determine 
Wo, the additional unknown dimensionful constant in tmChPT, if W is already determined 
elsewhere. 



^ In the tmChPT counting we use, these vertices are 0{p'^), and one needs at least two of them to contribute 
to the TTTT-scattering amphtude entering Liischer's formula given by Eq. (jlOp . Equivalently, one needs at 
least two such vertices to make a self-energy correction to the pion propagator that yields FVE^.. . 

^ The LEG W can be determined from the mass splitting between the charged and neutral pion [12]. But 
this is difficult in practice as calculating the neutral pion mass on the lattice involves quark disconnected 
contributions. 
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IV. NLO FORMULAE FOR THE PION MASSES AND DECAY CONSTANT IN 
FINITE VOLUME 

A. The standard contributions 



We provide here the complete resummed asymptotic formulae at NLO for the relative 
finite volume shift of the charged pion mass and decay constant. We split the contributions 

to the NLO FVEs for the charged pion mass into contributions in decreasing exponential 
importance: 

Rm± — Rm± (Aq) + Rm± (A±) , (16) 

where Rm± ( Aq) is the standard contribution due to neutral pions traveling around the whole 
volume (~ e"'^"), and Rm^{^±) that due to the charged pions (~ e~^^). 

The standard contributions Rm± start at 0{p'^), and it is easy to obtain the next order 
correction by evaluating TTTT-scattering at NLO in tmChPT.® Inserting this into the resummed 
Liischer formula, we have: 



Rm±{\) 



2A 



± 



E 



'm{n) 



(4) 



, (17) 



where i G {±, 0} is the isospin index, m{n) is the multiplicity of the integer vector n of 
length n= |n|, and 



7-(2) 

'm,o 



j(4) _ 

-'m,o ~ 



-B\ 



7-(2) 



0, 



(18) 



2h + - 
18. 

C OA 



5° + 
5° + 



20 8 
"9 
92 8 
¥ ~ 3 



B^ + -{Rl + 2Rl-m^o) , 



i?2 + -(lli?0-20i?^-32i?2) , 



with B^'' = B^'^i^Xi) and = i?J^'^(v^Ai) integrals given by 

— oo 

oo 

Ro ^'^ ^if^ Jdy Z/'e-V«[i+(2/M)^]A>^W (2 + 2iy) , for 



k even 
k odd 



(20) 



^ For simplicity, wc have set M^o = M^± in the NLO TTTT-scattering amphtude, as the effects of the charged- 
neutral pion mass splitting is higher order. The same goes for the amplitude entering the Liischer formula 
for F„. 
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where = MT^i/MT^± and 



(7-1 



^(x) = (7 log ^— + 2, a{x) = ^/l^^I/x . 

O" + i 



(21) 



For completeness we give here also the full NLO resummed asymptotic formula for FVEs 
in the neutral pion mass, although these will not be used in our numerical analysis. The 
formula reads: 

Rmq — Rmo{\) + Rmq{^±) , (22) 

where 

C± v-^ m{n) 



Rmo{K) 



2Ao^ 

n=l ^ 



(23) 



The integrals j can be expressed in terms of the I^^^^ integrals defined for the charged 
pion mass: 



^Molo(VnAo) = /i;_^±(v^Ao) -/]^,o(v^^o) 
= 2/£L(V^A±) . 



(24) 



For the pion decay constants, the decomposition of the finite volume shifts is similar to 
that given for the pion mass. We have: 



Rf± — -Rf±(Ao) + -Rf±(A±) 



(25) 



where 



and 



R-F± (Aj) 



n=l 



n 



41(V^A,)+e±/gi(V^A,) + 0(e|) 



(4)/ 



(26) 



7-{2) _ 7-{2) _ _ oO 
^F,0 ~ ^F,± ~ ^ 



7-{4) _ 



1 2- 

- + -£i - £4 

9 3 



5" + 



20 8 



8 4-4- 

— + + -£2 - 24 

9 3 3 



9 3 
5° + 



52 + ^ (2it:° + ml - %Rl - Rl' -2Rl' + ARl') , 
9 3 



+ 
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2i?° - 8i?^ - 32i?^ _ i?0' + ioi?i' + IQRl' 



(27) 



The function g{x) is related to the standard J one-loop function through g{x) = IGtt J{xM^±) 
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B. The parity-violating cubic interaction contributions 



As discussed above, because of the exponential enhancement, contributions due to parity- 
violating cubic interactions in tmChPT may have an effect at NLO in addition to the 
standard contributions to FVEs (arising from parity-conserving quartic interactions), despite 
being formally higher order. For the charged pion mass, they are given by 



2 ^± 



n=l * 



A J + 

— oo 



(28) 



Recall that w = M^o/(2M^±) < 1. As explained above, the leading exponential contri- 
butions to pion mass in finite volume arises first at 0{p^) in the chiral counting. Thus, 
contributions to the relative shift given here is 0{p^). Note that the sign of Ai?M± is fixed, 
and is opposite to the standard contributions given in Eq. ( 1T6|) . which tend to make the pion 
masses larger; ARm± would hence reduce this enlargement. 

For the pion decay constants, the exponentially enhanced contributions arise from the 
parity- violating interaction with the axial current. These have the same form as that for 
the charged pion mass except for a different pref actor involving the LEG Wiq. 

/\Rf± = — —ARm± , = . (29) 

Note that ARp^ is actually proportional to e, and unlike the case for pion mass above, its 
sign is not fixed because Wiq is unknown. We also remark that this contribution to the 
FVEs depends on how exactly is calculated. If the pseudoscalar density is used instead 
of the axial current as is the case here, ARp^ would be proportional to W instead of Wiq. 
This is a consequence of the fact that the Ward identity 

G^rrig = FMl (30) 

is violated by discretization effects in tmChPT. 



V. NUMERICAL ANALYSIS 



As a first illustration of the importance of the full NLO discretization effects in finite 
volume corrections, we evaluate them numerically using the formulae derived in the previous 
section for a box of size 2 fm and with pion masses close to those in the ETMC ensemble 
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Relative shift 




LO 




NLO 


Cubic interactions 


(%) 


ChPT 


tmChPT 


ChPT 


tmChPT 






0.38 


0.71 


0.74 


1.24 


-0.08 




-1.5 


-2.2 


-2.1 


-3.0 


0.06 



TABLE I: Relative finite volume corrections in percentage for a L = 2 f m box. The ChPT columns 
denote the case of degenerate pion with M^± = M^o = 0.33 GeV, while the tmChPT columns 
that where charged and neutral pions are split with M^± = 0.33 GeV and M^a = 0.27 GeV. At 
these mass values, F.,^ is about and so taken to be 0.11 GeV [^]. The 37r- vertex contributions are 
obtained by setting e = ^ and s5/e = 1. 



Bl 







The results are shown in Table HI where we give not only the LO and NLO standard 
contributions, but also the parity-violating cubic interaction contributions. As discussed 
above, these involve unknown LECs. So to evaluate them, we set e = ^ and s5/e = 1, which 
is a conservative choice. Crucially, we see when compared to the standard contributions, 
those arising from the parity-violating cubic interactions turn out to be subdominant after 
all, despite the exponential enhancement. Note as observed above, in the case of the 
two types of contributions have opposite signs so that finite volume corrections are reduced 
(albeit only slightly). In the case, the sign is chosen when we set s5/e = 1; it is otherwise 
free in general. 

Overall, Table |T] shows that the corrections to the continuum formulae are substantial, 
expecially for the pion masses, and this argues for a more detailed analysis of the ETMC 
data. Using the most recent ETMC data [8], we performed a new analysis to evaluate the 
impact on the finite volume corrections when effects due to the lighter neutral pion are 
included. A global fit to the ensembles Bi — B4, Bq, B7, Ci — C3 and Di is performed, 
which all have a charged pion mass below 500 MeV and a box size of at least 2 fm (so that 
ChPT can be safely applied).^ Our fit is purely statistical and does not include systematic 
errors, which have not been released. We recall that ensembles in the B (D) set have a 



Ensembles C5 and D2 have L ~ 1.6 fm, and thus too small a volume to have the FVEs reliably described 
by our formulae. Nevertheless, we have checked that including them in our fits does not deteriorate the 
quality of the fits, nor changes any of our conclusions. 
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coarser (finer) lattice than the C set. 

The forms we use to fit the ETMC data on the charged pion masses and decay constants 
at maximal twist read: 



M^±{a,L) = VM' 



F^{a,L) = F 



M' 



327r2F2 



1/2 



[l + R 



M±) , 



M' 



167r2F2 



[1 + Rf^)- 



(31) 
(32) 



Notice that we have not included ARm^ and ARp^ in our fitting forms: we have checked that 
the data at their current precision, have no sensitivity to these higher order parity-violating 
effects.^ We have however included parameters and Df in our fit, which account for the 
relative 0{a^) effects in and F^^ respectively. Although such effects are formally higher 
order (NNLO) in the counting we use, they may appear spuriously if the maximal twist 
condition, u = tt/2, is only determined to 0{a) accuracy. They are therefore included to 
provide useful diagnostics. 

In all, our fitting parameters at maximal twist thus consist of the tmChPT parameters 
(in lattice units) acF, 2acBQ, acA^, acA^ and a'^K, the ratios of lattice spacings aB/ac 
and a£,/aci and additional lattice discretization parameters a^D^ and o-cDf. The notation 
ax here denotes the lattice spacing for the ensemble set X, and we have chosen ac to be 
the reference lattice unit. The scale setting is accomplished by fixing = 92.4 MeV at the 
point where the ratio M-j^/Ft^ assumes its physical value. Note that the r^/ax data provide 
constraints on the ratios of lattice spacings. In our fitting the values of r^/ax are taken in 
the chiral limit as provided in Ref. 8|. We do not fit the LECs £1^2? which appear at NLO 
(or (9(p^)) in the finite volume corrections. They are fixed as in Ref. |l5| with their mass 
independent parts set to 



Tphys 



-0.4 ±0.6, 



iphys 



4.3 ±0.1. 



(33) 



We have investigated in our fitting, the effects of turning on the pion mass splitting [K) 
and higher order 0{a?) discretization effects {Dm and Dj) in various combinations. We 



^ When fitting with IS.Rm± and lS.Rp_^ included, we find that e is driven to be vanishingly small, as the 
data do not favour a reduction of FVEs in the charged pion masses that Ai?M± necessarily brings about. 
Because of this, fluctuations in the pion decay constant data tend to drive 5 into a run-away increase, as 
it has to compensate for the smallness of e (see Eq. ([211)), and this is seen in the fitting. As a result, no 
firm conclusions can be drawn about the presence of these parity- violating contributions from the data. 
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Fit 


I 


II 


III 


IV 


V* 


Pn 


O .tit: I X4: I 


O.Olt ID J 


O.^tJt -LO 1 


o.oo\±o J 


O.OOl ± -L I 


L 


4.69(3) 


4.63(4) 


4.75(5) 


4.70(6) 


4.71(6) 




1.0748(7) 


1.0736(9) 


1.0760(10) 


1.0749(11) 


1.0752(11) 




1 n28H'12'l 


1 n269('14') 


1 n283('12') 


1 0273n4'l 


1 0276l'10'l 


ar ffml 


0.0665(6) 


0.0676(7) 


0.0644(22) 


0.0649(28) 


0.0647(24) 


Or ffm) 

^ ID 'yi-LLi.j 


0.0856(6) 


0.0871(9) 


0.0811(40) 


0.0815(48) 


0.0813(43) 




05?8('5'> 


0536('6'l 


051 Q(^ 1) 


0523n 5"! 


0521 (^ A) 






0.0030(4) 




0.0032(11) 


0.0028(4) 








0.055(28) 


0.055(38) 


0.057(35) 








0.038(32) 


0.045(44) 


0.044(39) 




37.5/16 


32.6/15 


20.6/14 


16.1/13 


21.8/19 


p- value 


0.00 


0.01 


0.11 


0.24 


0.29 



E II: Results of all our fits of the ETMC data from ensembles Bi — B4^,Bq,Bi,Ci — C3 and 

1 Is]. 



TAB 
D 



The asterisk besides the fit number indicates the inclusion of the M^o data. The quantities 
M,r and F-,^ are physical pion mass and decay constant, and is the value of the quark mass 
corresponding to the physical pion mass. 

have also investigated the impact of including the M^o data on our fit. We do not include 
the MjrO data d priori because it is of a much lower quality compared to the M^± data (the 
uncertainty associated with it is at least an order of magnitude larger), and it has rather 
different systematics. Since the finite volume corrections are expected to be (9(1%), which 
would be easily subsumed by the C(10%) error in M^o, we do not apply them to the M^o 
data when including them in our fitting.^" Nevertheless, this should be done when the 
neutral pion mass can be calculated more reliably in the future, and we have provided the 
full NLO resummed formula for FVEs associated with M^o in Sec IIVI 

The results of our fits are shown in Table HTl To establish a baseline, we performed Fit I 
with pure ChPT fitting forms that included neither twisted mass nor NNLO discretization 



We have checked explicitly that applying finite volume corrections to the M^o data do not alter our fit in 
any way. 
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effects, and we see it gives the worst description of the data. Turning on just the pion mass 
sphtting in Fit II induced only a marginal improvement, as the decreased by 2.4% with 
the number of degrees of freedom, naof , reduced by one. However, turning on instead and 
Df in Fit III produced a dramatic improvement with a 45% reduction in the compared to 
Fit I as ndof reduced by two. A further improvement is gained when all parameters associated 
with discretization effects, K, Dm and Df, are simultaneously turned on in Fit IV, as the 
decreased another 22% compared to Fit III with ridof reduced by one. Finally the best fit is 
Fit V* when M^o data are also included (indicated by the asterisk beside the fit number). 

For each fit, we have also calculated the p- value, namely the probability of obtaining a 
normalized greater than that actually found from the fit. It is particularly illuminating 
to compare the p- values in Table [TTl as it shows that despite the substantial improvement in 
X^, Fit III is not yet fully convincing in terms of its p- value. Only after including the twisted 
mass effects as in Fit IV does the p-value rise to a much more acceptable level. Finally, the 
inclusion of the M^o data gives a further increase in the p- value in Fit V*. 

We remark here that for the pion mass-squared splitting, a^K = M^± — M^o = AM^, 
its fitted value changes little whether Mt^o data are included in the fit or not. Incidentally, 
including M^^o data seems to reduce the uncertainty in a^K substantially, and leads to a 
result which is (statistically) different from zero by about 7cr compared to a bit under 3a 
when M^o data are not included. By combining the analysis of the FVEs of the charged 
pion masses and the M^^o data, one can see unambiguously the mass splitting predicted 
in tmChPT, and that AM^ > 0. This represents the first determination of the tmChPT 
quantity WqW from tmLQCD simulations. We emphasize here that AM^, and hence M^o, 
can already be determined from the FVEs alone. Although large uncertainties are associated 
with the determination of M^o either through analyzing FVEs or by direct lattice calculation, 
it is reassuring to see that both provide compatible results. 

To provide further insight into exactly where improvements arise by taking into account 
twisted mass effects, we give a breakdown of the individual contributions to the total 
from each ensemble for Fit III and IV in Table lllli We would like to see the individual 
contribution from each and F^^ data, and this is given by x%i^ and x'f^ in Table lllli Note 
that since and F^^ are correlated non-trivially in any given ensemble, Xm^ and Xf^ do 
not sum up to the full contribution from that particular ensemble, Xpair? which takes into 
account the correlation between the pair: they are x^ values calculated for each and F„ 
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Ensemble data 






Fit III 






Fit IV 




Ensemble 


L/a 


aflq 




xk 


X^ • 




xk 


Apair 


i>l 


A 

z4 


0.0040 


3.16 


0.51 


3.24 


0.63 


0.42 


2.23 


B2 


24 


0.0064 


1.84 


0.00 


3.74 


1.91 


0.00 


3.68 


TD 


OA 


0.0085 


0.03 


0.39 


0.40 


0.00 


0.20 


0.23 


TD 

J=>4 


z4 


0.0100 


1.61 


0.40 


1.87 


1.52 


0.37 


1.76 


7-) 


62 


0.0040 


1.37 


0.53 


1.37 


0.19 


0.03 


0.50 


B7 


32 


0.0030 


0.02 


0.57 


1.64 


0.74 


0.01 


1.87 


Ci 


32 


0.003 


4.85 


1.04 


4.88 


3.76 


0.61 


3.76 


C2 


32 


0.006 


0.16 


0.06 


0.38 


0.18 


0.02 


0.30 


C3 


32 


0.008 


0.00 


1.25 


1.46 


0.07 


0.66 


0.99 


Di 


48 


0.0020 


1.20 


0.96 


1.51 


0.58 


0.49 


0.75 



TABLE III: The contribution of each individual ensemble to the total for Fit III and IV. The 

quantities and xk denote the contribution each and i<V data makes separately to the 

total x^ neglecting the correlation between them. The full contribution from the correlated pair is 
denoted by Xpair- 

data separately neglecting the correlation 

A comparison of Fit III and IV shows very clearly for which ensemble and for which 
quantity is the improvement from taking into account twisted mass effects the most signifi- 
cant. The largest improvements are seen in M^r for ensembles Bi and Ci, which indeed have 
small volumes and small pion masses. For these two ensembles, the relative improvement in 

is also quite evident. In all, we see clear improvements across-the-board. 

VI. CONCLUSIONS 

In this paper we have analyzed the finite volume corrections to pion masses and decay 
constant in tmChPT. The presence of a lighter neutral pion gives exponentially enhanced 

finite volume corrections, and we have calculated these to NLO accuracy. Wc have performed 
a detailed analysis of the ETMC data and found out that they are better fitted by the 
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formulae derived here than by the continuum ones j4|, particularly in regards to the finite 
volume dependence. We are able to extract the pion mass splitting predicted in tmChPT 
from analyzing FVEs on charged pions alone without having to calculate directly the neutral 
pion mass. An important benefit of this is, as far as we know, a first determination of the 
LEG W of tmChPT. This example shows that, though small, FVEs can be successfully 
used to determine interesting physical observables. 

Other LECs of tmChPT appear in parity-violating cubic interactions, which give expo- 
nentially enhanced FVEs. Despite the fact that they are formally of higher order in the 
combined twisted mass chiral expansion, we have calculated their analytical form and in- 
vestigated their effects in our fitting in the hope that the exponential enhancement may be 
large enough to allow them to be seen, and thus enable extraction of more LECs of tmChPT. 
Unfortunately, their contribution is found to be small, and below the precision of the present 
ETMC data. 

As the precision of the lattice calculations increases and as the simulations move towards 
lighter pions, the corrections discussed here will become even more important and may have 
an impact also on the extracted physically relevant parameters. We suggest that future 
tmLQCD lattice studies perform analyses of FVEs taking full account of the twisted mass 
discretization effects, as was done in deriving our formulae in this paper. 
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